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We examine the phase transition of polymer adsorption as well as the underlying kinetics of 
polymer binding from dilute solutions on a structureless solid surface. The emphasis is put on the 
properties of regular multiblock copolymers, characterized by block size M and total length A'^ as 
well as on random copolymers with quenched composition p of sticky and neutral segments. The 
I macromolecules are modeled as coarse-grained bead-spring chains subject to a short-ranged surface 

■ adhesive potential. Phase diagrams, showing the variation of the critical threshold for single chain 

' adsorption in terms of M and p are derived from scaling considerations in agreement with results 

' from computer experiment. 

Using both scaling analysis and numerical data from solving a system of coupled Master equations, 
we demonstrate that the phase behavior at criticality, and the adsorption kinetics may be adequately 
predicted and understood, in agreement with the results of extensive Monte Carlo simulations. 
Derived analytic expressions for the mean fraction of adsorbed segments as well as for Probability 
Distribution Functions of the various structural building blocks (i.e., trains, loops, tails) at time 
t during the chain attachment process are in good agreement with our numeric experiments and 
provide insight into the mechanism of polymer adsorption. 
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I. INTRODUCTION 



' ^ ' The adsorption of polymers on solid surfaces is a long-standing problem which plays an important role in a number 
of applications in technology (protective coatings, adhesives, lubricants, stabilization of colloids, flocculation, etc.) and 
biology (adsorption of biopolymers, etc.). As a phenomenon it poses a number of challengiiig scientific problems^, 
H, SB too. Important theoretical contributions have been made by BirshteinQ, de Gennesfq, Eisenriegler et al.f^. 
Later studies have covered adsorption of polyelectrolytes dynamics of adsorbed chains 9] and adsorption on 
chemicall y h eterogeneous surfaces [l0| . The close relationship between analytic theory and computer experiments in 
this field [lOl . [Til , [l^ . [Tsl . [T^ has proved especially fruitful and instructive. 
J — ] While the investigations mentioned above have been devoted exclusively to homopolymers, the adsorption of copoly- 
mers (e.g., random or multi-blocks copolymers) still poses open questions. Thus, for instance, the critical adsorption 
' potential (CAP) dependence on block size M at fixed concentration of the sticking j4-monomers is still unknown as are 
, the scaling properties of regular multi-block copolymers in the vicinity of the CAP. From the theoretical perspective, 
' the case of diblock copolymers has been studied by means of the Grand Canonical Ensemble (GCE)-approach [15l.[l^. 
0^ within the Self-Consistent Field (SCF)-approach t)y Monte Carlo computer simulations [l^, [l^]. The case 

of random copolymers adsorption has gained comparatively more attention by researcher so far. It has been inves- 
^ tigated by Whittington et al. [2l|, [2^ using both the annealed and quenched models of randomness. The influence 
of sequence correlations on the adsorption of random copolymers has been treated by means of the variational and 



replica method approach 23j . Sumithra and Baumgaertner [24| examined the question of how the critical behavior 
of random copolymers differs from that of homopolymers. Thus, among a number of important conclusions, the 
results of Monte Carlo simulations demonstrated that the so called adsorption (or, crossover) exponent 4> (see below) 
is independent of the fraction of attractive monomers n. 

The adsorption kinetics of polymers has been intensively studied both experimentally [1^, [l^ and theoretically [l^, 
[28l . [2^ . [sol , [sil . [3^ . [33I since more than two decades now. A key parameter thereby is the height of the free energy 
adsorption barrier that the polymer chain has to overcome so as to bind to the surface. High barriers are usually 
referred to as cases of chemisorption as opposed to those of physisorption which are characterized by low barriers for 
adsorption. Depending on the strength of the binding interaction e, one distinguishes then between weak physisorption 
when e is of the order of the thermal energy ksT (with ks being the Boltzmann constant), and strong physisorption 
when e > 2kBT. One of the important questions concerns the scaling of the adsorption time Tads with the length 
of the polymer chain N in dilute solutions. For homopolymers in the reg ime of strong physisorption (that is, for 
sticking energy considerably above the CAP) computer experiments [H, ISll ls^ suggest Tads oc N" where a is related 
to the Flory exponent va,sa = l + iyK 1.59. This result follows from the assumed zipping mechanism in the absence 
of a significant barrier whereby the chain adsorbs predominantly by means of sequential, consecutive attachment of 
monomers, a process that quickly erases existing loops on the substrate. For the case of weak adsorption, one should 
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mention a recent study [33|, where one finds in contrast a = (1 + 2i/)/(l + r/) « 1.37 which suggests shorter time 
scale for surface attachment. In chemisorption, the high barrier which attaching monomers encounter slows down the 
binding to the surface, the chain gains more time to attain equilibrium conformation and the adsorption process is 
believed to involve large loop formation giving rise to accelerated zipping mechanism [29l [soj . The predicted value of 
a in agreement with MC results is a ~ 0.8 ± 0.02 [s^ . A comprehensive overview of experimental work and theoretic 
considerations may be found in the recent review of O'Shaughnessy and Vavylonis [30l |. 

In the present contribution we focus on copolymer physisorption, extending thus the aforementioned studies of 
homopolymers statics and kinetics. We show how scaling analysis as well as different MC-simulation methods help 
understand the critical behavior of multi-block and random copolymers. It turns out that the critical behaviour of 
these two types of copolymers can be reduced to the behavior of an effective homopolymer chain with "renormalized" 
segments. For multi-block copolymers one can thus explain how the adsorption threshold depends on the block length 
M and even derive an adsorption phase diagram in terms of CAP against M. In the case of random copolymers, 
the sequence of sticky and neutral (as regards the solid substrate) monomers within a particular chain is usually 
fixed which exemplifies a system with quenched randomness. Nevertheless, close to criticality the chain is still rather 
mobile, so that the sequence dependence is effectively averaged over the time of the experiment and the problem can 
be reduced to the easier case of annealed randomness. We show that the MC-findings close to criticality could be 
perfectly described within the annealed randomness model. 

For both regular multiblock as well as for random copolymers, we compare the predicted kinetics of adsorption in 
the regime of strong physisorption, to consistent numeric data derived from simulations and coupled Master equations. 
We demonstrate that the observed adsorption kinetics is close to that of homopolymers and suggest interpretation of 
typical deviations. Eventually, we should like to stress that the complex polymer hydrodynamics near an interface 
has remained beyond the scope of this paper. 



II. SIMULATION METHODS 

Apart from the frequently used Bond-Fluctuation Method (BFM) [H, [s^l, two coarse-grained models, a bead- 
spring off-lattice model P and a cubic lattice model implementing the so called pruned-enriched Rosenbluth method 
(PERM) [31, are used to test theoretical predictions. 



A. Off-lattice bead-spring model 

In our computer simulations we use a coarse grained off-lattice bead spring model [ll| to describe polymer chains. 
The system consists of a single chain tethered at one end to a flat structureless surface so as to avoid problems with 
translational entropy depending on the box size. There are two kinds of monomers: "A" and "B", of which only 
the "A" type feels an attraction to the surface. The surface interaction of the "A" type monomers is described by a 
short-range square well potential Uw{z) = e ior z < 6 and Um{z) = otherwise, whereby the range 6 = 0.125 (in units 
of the maximal bond length extension Imax between adjacent beads). The effective bonded interaction is described 
by the FENE (finitely extensible nonlinear elastic) potential. 



Ufene = -K{1 - lofln 



l-lo 



(2.1) 



with K = 20, Imax = 1, ^0 = 0.7, l,nin = 0.4 

The nonbonded interactions are described by the Morse potential. 



exp(-2a(r - rmm)) - 2exp(-a(r - r™„)) 



(2.2) 



with a = 24, r™„ = 0.8, eM/ksT = 1. 

We use periodic boundary conditions in the x — y directions and impenetrable walls in the z direction. We study 
homopolymer chains, regular multiblock copolymers, and random copolymers (with a fraction of attractive monomers, 
p = 0.25, 0.5, 0.75) of length 32, 64, 128, 256 and 512. The size of the box is 64 x 64 x 64 in aU cases except for the 
512 chains where we use a larger box size of 128 x 128 x 128. The standard Metropolis algorithm is employed to govern 
the moves with self avoidance automatically incorporated in the potentials. In each Monte Carlo update, a monomer 
is chosen at random and a random displacement attempted with Ax, Ay, Az chosen uniformly from the interval 
—0.5 < Ax, Ay, Az < 0.5. The transition probabihty for the attempted move is calculated from the change AU of 
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the potential energies before and after the move a,s W = exp{—AU/kBT). As for a standard Metropohs algorithm, 
the attempted move is accepted if W exceeds a random number uniformly distributed in the interval [0, 1). 

Typically, the polymer chains are originally equilibrated in the MC method for a period of about 10^ MCS (depend- 
ing on degree of adsorption e and chain length N this period is varied) whereupon one performs 200 measurement 
runs, each of length 8 x 10® MCS. In the case of random copolymers, for a given composition, i.e., percentage p of the 
A— monomers, we create a new polymer chain in the beginning of the simulation run by means of a randomly chosen 
sequence of segments. This chain is then sampled during the course of the run, and replaced by a new sequence in 
the beginning of the next run. 



FIG. 1: Schematic representation of a grafted chain close to criticality. Snapshot of a regular multiblock copolymer with length 
iV = 2048 and block size M = 8 from the PERM simulation. 



The adsorption of a diblock AB copolymer with one end (monomer A) grafted to a flat impenetrable surface and 
with only the A monomers attractive to the surface is described by self-avoiding walks (SAW) of — 1 steps on a 
simple cubic lattice with restriction z > 0. The partition sum may be written as 



where Aff{Ns) is the number of configurations of SAWs with N steps having Ng sites on the wall, and q = e'^/'^^^ is 
the Boltzmann factor (e > is the attractive energy between the monomer A and the wall). Clearly, any copolymer 
will collapse onto the wall, if q becomes sufficiently large. Therefore one expects a phase transition from a grafted 
but otherwise detached phase into an adsorbed phase, similar to the transition observed for homopolymers. 

The pruned-enriched Rosenbluth method (PERM) p3 |. also used in our simulations, is a biased chain growth 
algorithm with resampling (" population control" ) and depth-first implementation. Polymer chains are built like 
random walks by adding one monomer at each step. Thus the total weight of a configuration for a polymer consisting 
of A^ monomers is a product of those weight gains at each step, i.e. Wn — ^it^^Wi. As in any such algorithm, there 
is a wide range of possible distributions of sampling so we have the freedom to give a bias at each step while the 
chain grows, and the bias is corrected by means of giving a weight to each sample configuration, namely, Wi Wi/pi 
where pi is the probability for putting the monomer at step i. In order to suppress the fluctuations of weights as the 
chain is growing, the population control is done by "pruning" configurations with too low weight and "enriching" the 
sample with copies of high-weight configurations. Therefore, two thresholds are introduced here, = c^Zn and 
W~ = Zn, where Z„ = -p- '^^.^njig Wn from the Af„ trail configuration is the current estimate of partition sum 
at the n — 1 step, c"*" and c~ are constants of order unity and c"'"/c~ ~ 10. In order to compare with the results 
obtained by the first MC method, we simulate homopolymers of length TV = 2048 and multi-block copolymers with 
block size A/ = 2*^^ , fc = 0, 1,2, • • • ,9 (see Fig. [T]). The number of monomers is increased to A^ = 8192 as the 
block size increases. There are 10^ 10® independent configurations for each measurement. We also simulate random 
copolymers of A^ = 2048 monomers with composition p = 0.125, 0.25, 0.50, and 0.75. 




B. Coarse-grained lattice model with PERM 




(2.3) 
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III. SCALING BEHAVIOR AT CRITICALITY 



A. A homopolymer chain 

It is well known 0, [13, that a single polymer chain undergoes a transition from a non-bound into an adsorbed 
state when the adsorption energy e per monomer increases beyond a critical value Cc ~ ksT (where T stands for the 
temperature of the system) . The adsorption transition can be interpreted as a second-order phase transition at the 
critical point (CAP) of adsorption e = ec in the thermodynamical limit, i.e., N ^ oo. Close to the CAP the number 
of surface contacts Ng scales as Ns{e = ec) ~ N'^. The numerical value of (j) is somewhat controversial and lies in a 
range between = 0.59 (ref. ff\) and (j) = 0.484 (ref. [lj|), we adopt however the value = 0.50 ± 0.02 which has 
been suggested as the most satisfactory [l^ by comparison with comprehensive simulation results. 

How does polymer structure vary with adsorption strength? Consider a chain tethered to the surface at the one 
end. The fraction of monomers on the surface n — N^/N may be viewed as an order parameter measuring the degree 
of adsorption. In the thermodynamic limit N — > oo, the fraction n goes to zero (« 0{1/N)) for e ^ ec, then near 
Cc, n ^ N'^~^ whereas for e ^ ec (in the strong coupling limit) n is independent of N. Let us measure the distance 
from the CAP by the diniensionless quantity k = (e — ec)/ec and also introduce the scaling variable rj = kN'^. The 
corresponding scaling ansatz^S^] is then 71(77) = N"^^^ G {rj) with the scaling function 

G{rj) = h'''\^, ' ^""^ ^^'^ (3.1) 
^ I , for 77 > 1 ^ ' 



The resulting scaling behavior of n follows as. 



n (X < 



l/N , for K<0 

7V^-i , for K^O (3.2) 

Kfi-"^)/"^ , for K > 1 



The gyration radius in direction perpendicular to the surface, Rg±(ri), has the form Rg±(j]) ~ aN'^Qg± (77). One 
may determine the form of the scaling function ^?g_L(77) from the fact that for k < one has Rg± ^ aN'^ so that 
Qg±^ = const. In the opposite limit, 77 ^ the A^-dependence drops out and Gg^iv) ~ 77^"^/"^. In result 

^^^(^)°^\."''/^ , for 7;»0. ^'-'^ 

The gyration radius in direction parallel to the surface has similar scaling representation, i?g||(77) = aN'^Qg\\ {rj). 
Again at k < the gyration radius ~ aN" and Gg\\ — const. At 77 ^ the chain behaves as a two-dimensional 
self-avoiding walk (SAW), i.e., ~ aN'^^, where 1^2 — 3/4 denotes the Flory exponent in two dimensions. In result, 
the scaling function behaves as Ggwirj) = 77'"^^'^^/''^, at 7/^0. Thus 

fi. f ^ /°^" , at 77<0 

The study of the ratio r(77) = Rg^/Rg\\ of gyration radius components is a convenient way to find the value of ec 
In fact, from the previous scaling equations r{ri) = Gg±{ri)/Gg\\{v)- Hence at the CAP, i.e., at 7; ^ the 
ratio r(0) = const, is independent of N. Thus, by plotting r vs. e for different N all such curves should intersect at a 
single point which gives ec - cf. Fig. [2k. 

If the well-known picture of blobs [6l| is invoked, then in the limit KiV^ ^ 1 the adsorbed chain can be visualized 
as a string of adsorption blobs which forms a pancake-like quasi-two-dimensional layer on the surface. The blobs are 
defined to contain as many monomers g as necessary to be on the verge of being adsorbed and therefore carry an 
adsorption energy of the order oiksT each. The thickness of the pancake Rg± corresponds to the size of the blob while 

the chain conformation within a blob stays unperturbed (i.e., it is simply a SAW), therefore g ~ {Rg±l a)^^^ = n^^/'^ 
where we have used Eq. 13.31 The gyration radius can be represented thus as 

Rg\\ = Rgi_ ' oc K^^^-'^/'^TV''^ (3.5) 

and one goes back to Eq. 13.41 which proves the consistency of the adsorption blob picture. Generally speaking, the 
number of blobs, N/ g ~ hi^/'^N , is essential for the main scaling argument in the above-mentioned scaling functions. 
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FIG. 2: Variation of the ratio Rg±/Rg^^ for a copolymer of block size Af = 2 (a), and of the specific heat per monomer, Cv{t)/N 
for a homopolymer (b) with surface energy e for diflPerent chain length N. One gets ef =^ « 0.556 and « 0.284 (PERM). 



The adsorption on a plane at k > is due to free energy gain which is then proportional to the number of blobs, i.e., 
F — i^buik c< ^N/ g ~ ^Hi^/'^N . The expression for the specific heat per monomer may be then obtained as 

9^(i^-^bulk) _a f„ 

Cv = (X K " , (3.6) 

where a = 2 — (t>~^ ■ If </> = 0.5 then a = and the specific heat does not diverge but rather undergoes a jump at the 
CAP - cf. Fig. [13. 



B. Multiblock Copolymer Adsorption 

One may now consider the adsorption of a regular multi-block copolymer comprising monomers A which attract 
(stick) to the substrate and monomers B which are neutral to the substrate. In order to treat the adsorption of a 
regular multi-block ^i3-copolymer it may be reduced to that of a homopolymer which has been considered above. 
Thus a regular multi-block copolymer can be treated as a "homopolymer" where a single Ai3-diblock plays the role 
of an effective monomer [s^ . Let each individual diblock consist of an attractive A-block of length M and a neutral 
B-block of the same length. Upon adsorption, the A-block would form a string of blobs whereas the i?-part forms a 
non-adsorbed loop or a tail. The free energy gain of the attractive block may be written then (in units of ksT) as 
Fattr = -^iZ-^M where k ee (e - ei')/ej now measures the normalized distance from the CAP ej? of a homopolymer. 
The neutral i3-part which is most frequently a loop connecting adjacent A-blocks, but could also be a tail with the 
one end free, contributes only to the entropy loss i^rop = (7 — 7ii) InM where the universal exponents 7 and 711 are 
well known pjq (e.g. in 3_D - space 7 — 1.159, 711 = —0.390). If a tail is involved, one should also use the exponent 
71 — 0.679 albeit this does not change qualitatively the expression for i^rcp- These expressions reflect the standard 
partition functions for a free chain, a chain with both ends fixed at a two points, and for a chain, tethered by the one 
end [3^ . In result the effective adsorption energy of a diblock is [l^l : 

E{M) = K^I't'M - (7 - 711) InM . (3.7) 

One can tackle the problem of regular copolymer adsorption by mapping it on that of a "homopolymer" , consisting of 

TV = N /2M such effective units by using a — > aM^ , k — > A = "^~f , N — > A/" where a denotes the monomer 

size, and E'^ is the critical adsoption energy of the renormalized homopolymer. Generally, one would expect E'^ to 
be of the order of e^, reflecting the model dependence of the latter. At the CAP of the multiblock chain one has 
A = 0, thus one can estimate the deviation k^^ , of the corresponding critical energy of adsorption, e*^, from that of 
a homopolymer, namely 

,,M- ^f -4 / (7-7ii)lnM + gn ^ 

^c=^^=( j-j ■ (3.8) 
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FIG. 3: (a) Phase diagram showing k*^ = [e^' — tc)/^c plotted vs 1/M for multiblock copolymers with various values of M 
and = 1.716. The curve gives the best fit of Eq. 13.81 (b) The CAP, e^, plotted vs the composition p for random copolymers. 
Symbols denote the CAP for multiblock copolymers with block size M. The cureve corresponds to the best fit of Eq. 13.131 
Data obtained employing the ofi'-lattice bead-spring model. 



where we have used Eg. 13.71 It is seen from the phase diagram, Fig.[3K, that the deviation k^^ , Eg. 13.81 with (p = 0.5, 
steadily grows with decreasing block length M, in agreement with the computer experiment. The fraction of effective 
units on the surface obeys the same scaling law as given by eg [312 i-e., n = ^ = N't'-^G (AAf^) which becomes 
accurate, provided (i) k <C 1 but M '3> I such that InAf ^ 1 and k^^'^M ^ 1, and (ii) TV ^ 1. Thus, within each 
effective unit of the A-monomers only Ms will be adsorbed at criticality whereby scales as Ms = M'^G (kM"^) 
so that the total number of adsorbed monomer Ng = NsMg = NsM'^G {^kM'^^. The adsorbed fraction of monomers 
then is expected to scale with both N and M as 

n cx N^-^ G {kM^) G [ A ^ J . (3.9) 



For sufficiently strong adsoption, k\/M ^ 1 and A a/ N/AI ^ 1, one gets thus n oc kA 

The gyration radius compon 
TZ± ~ aA-'^/'f'M'', which yields 



The gyration radius component in direction perpendicular to the surface 'R-g± = aN"^ Qg± f A(^)'^j becomes 



aM-'E^ 

n± ~ £ 5- 3.10 

[k2M- (7-711) In M-S^f" 



In a similar manner, the gyration radius component parallel to the surface has the form T^^y = aN^ Qg\\ (^)'' 
which in the hmit A^N/M > 1 results in 7^g|| ~ a ( ) A^''^ i.e.. 



a[K2M-(7-7ii)lnM-i?,^]'^"^~'^^ ^ , ^ 

7^„|| ~ — ^ 'i^ "-^ N""^ . 3.11 

As shown in Fig. [U one finds indeed the expected scaling behavior which is demonstrated by the collapse of the 
simulation data on few "master curves" , absorbing cases of different strength of adsorption k. Like in the homopolymer 
case, one can define a blob length Q ^ {TZ^_/af'^" ^ A^^/'^M which in the strong adsorption limit, A > 1, approaches 
the block length, Q ^ M , as it should be. 



C. Random Copolymer Adsorption 

The adsorption of a random copolymer on a homogeneous surface has been studied by Whittington et al. [2l|, [l^l 
within the framework of the annealed disorder approximation. Physically this means that during the measurements 
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FIG. 4: Log-log plots of the order parameter n (a), and iigy/A^^" and Rgl^/N'^" vs kTV"^ (b) with u = 0.588 and U2 = 3/4. 
The straight lines indicate the asymptotic behaviour of the scaling functions given by Eqs. (|3.9p . (|3.10p . and 13. Ill for regular 
multi-block copolymers with block size M — 16 (PERM). 



the chain touches the substrate at random in such a way that one samples ah possible monomers sequences along the 
backbone of the macromolecule. Following this assumption [2l[, let CAr(n) be the number of polymer configurations 
such that n units have contact with the surface simultaneously. The percentage of ^-monomers (composition) is 
denoted by p. In the annealed approximation one then averages the partition function over the disorder distribution, 
i.e., 



N 



= E E p"-(l-p)"-"''e™'' 

N N 



(3.12) 



where e^g is the attraction energy of an effective homopolymer. From Eq. l3.12l one can see that the annealed problem 
is reduced to that of a homopolymer where the effective attractive energy is defined as e^g = In [pe'^ + 1 — p]. Since the 
homopolymer attraction energy at the CAP is e'^g = ejf, the critical attraction energy e = of a random copolymer 
will be 



In 



exp ej? -f p — 1 



P 



(3.13) 



where the composition < p < 1- Atp— >0 oo whereas at p — 1 e^? = ej?. This prediction which has been 



recently confirmed by Monte Carlo simulations [37[, is plotted in Fig. [3Jd. It shows that close to criticality the chain 
is still rather mobile, so that the sequence dependence is effectively averaged over the time of the experiment and the 
of quenched disorder can be reduced to that of annealed randomness. 



IV. ADSORPTION KINETICS 



A. Variation of the adsorbed fraction with elapsed time - Theory 

We illustrate here how one can use the "stem - flower" notion of adsorbing linear macromolecule, suggested by Descas 



et al. [3l|, to describe the observed "zipping" dynamics [33] of adsorption not only in terms of the average fraction of 
adsorbed segments but to include also time-dependent train and tail distribution functions as main constituents of 
the dynamic adsorption theory. The simple "zipping" mechanism along with the underlying stem-flower model are 
illustrated in Fig.[S^,b. The number of the adsorbed monomers at time t is denoted by n(t). The nonadsorbed fraction 
of the chain is subdivided into two parts: a stretched part ("stem") of length m{t), and a remaining part ("flower") 
which is yet not affected by the tensile force of the substrate. The tensile force propagation front is at distance R{t) 
from the surface. The rate of adsorbtion is denoted as v{t) = a '^^J'P , where a is the chain segment length. 



8 



A single adsorption event occurs with energy gain e and entropy loss In^fi^ / (12) , where ^3 and fj,2 are the connectivity 
constants in three and two dimensions, respectively [s^, so that the driving free energy F^^. — e — ksTlndi^/ 112) 
whereas the driving force /dr = F^j-fa. The friction force is related to the pulling of the stem at rate v{t), i.e., 
/fr = Co a m{t) ^^^^ where Co is the Stokes friction coefficient of a single bead. The equation of motion, following 
from the balance of driving, and drag force, is then /^r = Inspection of Fig. [5^ suggests that the distance R{t) 
between the flower and the plane changes during the adsorption process until the flower is eventually "consumed" . In 
so doing R{t) obeys two relationships: R{t) « a[n{t) + m{t)]'^ (because it is actually the size which the chain portion 
n{t) + m{t) occupied before the adsorption has started) and R{t) m m{t) (up to a geometric factor of order unity). 
Thus n{t) w ■m{ty^'^ — m(<) which yields m(t) « n{tY for the typically long stems ■m{t) ^ 1. From the resulting 
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Monomer index 



FIG. 5: (a) Chain conformation at successive time moments during the adsorption process for a polymer with A'' = 256. The 
z-coordinate of the i-th monomer is plotted against monomer index i. (b) Stem-flower picture of the adsorption dynamics. The 
total number of adsorbed monomers at time t is denoted by n{t). The tail which, contains all nonadsorbed monomers, consists 
of a stretched part, a "stem", of length m{t), and of a nonperturbed part - a "flower". The rate of adsorption is v{t). The 
distance between the surface and the front of the tension propagation is R{t). 

equation Co njtY dnjt) / dt = fdn/o? then follows n{t) oc « t^-^^ which is in good agreement with simulation 

results [2i,[3ll|3i]. 



B. Time evolution of the distribution functions - Theory 



Consider the instantaneous number of adsorbed monomers n at time t, i.e., the total train length) distribution 
function P(rt, t). Using the 'Master Equation' method [ssj . one may derive a system of coupled kinetic equations for 
P(n, t) by treating the zipping dynamics as a one step adsorption / desorption process within an elementary time 
interval. Assuming that the corresponding rate constants ?«^(n), 'w~{n) of monomer attachment / detachment are 
related by the detailed balance condition [38| (which is an approximation for a non-equilibrium process), one can 
flx their ratio u;+(n — \)/w^[n) = exp(_Fdr/fcBr), and even fully specify them by introducing a friction-dependent 
transmission coefficient q{Tn) = ksT/ (a^C,) = ksT / {a^QQm) (whereby the stem length m depends on the total train 
length n, according to n « w}/^ — m). Then the one-step Master Equation reads [3^| 

— P{n,t) = w-(n+l) P(n + l,t) + w+{n-l) P(n-l,t) 
dt 

~ w+{n) P{n,t) -w^{ji) P{n,t) , (4.1) 
which along with the boundary conditions 

j^P{l,t) = w-{2)P{2,t)-w+{l)P{l,t), forn-1 (4.2) 

^P{N,t) = w+{N-l)P{N-l,t)-w-{N)P{N,t), forn = iV 

and P{n, t — 0) — 6{n — I) fully describe the single chain adsorption kinetics. 
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The equation of motion for the mean number of adsorbed segments {n) 
Eq. ((i7T|) . assuming for simphcity P{N,t) = P{0,t) = 0: 



^,^-1^ nP(ri, t), can be obtained from 



- {n) = -{w {n)) + {w+in)) 
With the relations for the rate constants, w^{n), w^^n), this equation of motion becomes 



Co — n{t) 



knT r 



-Far/kBT 



(4.3) 



(4.4) 



where for brevity we use the notations n(t) — (n) and m(t) ~ (m). Note that Eq. (|4.4p reduces to the kinetic 
equation 31], derived at the end of Section flV Al for weak driving force, F^r ^ ksT, by neglecting fluctuations in 
the zipping mechanism. Evidently, by taking fluctuations into account, Fdr/a, is replaced by a kind of effective second 
virial coefficient {kBT/a)[l — exp{~Fdr/kBT)]. Thus, the zipping as a strongly non-equilibrium process cannot be 
treated quasistatically by making use a simple "force balance" . 



C. Order parameter adsorption kinetics - MC results 

The time variation of the order parameter n(t)/N (the fraction of adsorbed segments) for homopolymer chains of 
different length N and strong adhesion e/ksT = 4.0 is shown in Fig. [6^,b whereby the observed straight lines in 
double-log coordinates suggest that the time evolution of the adsorption process is governed by a power law. As the 
chain length N is increased, the slope of the curves grows steadily, and for length N = 256 it is equal to k, 0.56. 
This value is close to the theoretically expected slope of (1 -I- v)~^ ~ 0.62. The total time r it takes a polymer chain 
to be fully adsorbed is found to scale with chain length as r oc iV" whereby the observed power a^, 1.51 is again 
somewhat smaller than the expected one 1 + v k, 1.59, most probably due to finite-size effects. One may also verify 
from Fig. [BJd that for a given length N the final (equilibrium) values of the transients at late times t — *■ oo grow while 
the curves are horizontally shifted to shorter times as the surface potential gets stronger. Nontheless, the slope of the 
n{t) curves remains unchanged when e/ksT is varied, suggesting that the kinetics of the process is well described by 
the assumed zipping mechanism. The changing plateau height may readily be understood as reflecting the correction 




t [MCS] t[MCS] 

FIG. 6: (a) Time evolution of the order parameter (fraction of adsorbed segments) for four different chain lengths TV = 
32, 64, 128, and 256 at surface potential e/ksT = 4.0. The slope of the = 256-curve is 0.56. The inset shows the scaling of 
the adsorption time with chain length, r oc A^^'^^. The time r is determined from the intersection point of the late time plateau 
with the tangent to the respective n(t)-curve. (b) Adsorption kinetics for different strengths e of the surface potential. 
The variation of the plateau height (i.e., the fraction of adsorbed monomers at equilibrium) with e is depicted in the upper 
inset where the solid line nt^oo = 1 — 5exp {—e/ksT) describes the equilibrium number of defects (vacancies). The lower inset 
shows a collapse of the adsorption transients on a single 'master curve', if the time axis is rescaled appropriately. 

in the equilibrium fraction of adsorbed monomers due to the presence of defects (vacancies) for any given value of 
e/ksT. For the transients which collapse on a master curve, cf. the second inset in Fig.O), one may view the rescaling 
of the time axis by the expression t — > i[l — 13.7 exp -jr^] as a direct confirmation of Eq. 14.41 where the time variable 
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t may be rescaled with the driving force of the process (i.e., with the expression in square brackets). The factor 
« 13.7 gives then the ratio fJ-3/fJ.2 of the effective coordination numbers in 3- and 2-dimensions of a polymer chain 
with excluded volume interactions, /is and fi2 are model-dependent and characterize, therefore, our off-lattice model. 




, , 1 1 lU 3 4 , , , 

lo' lo" lo' 10 10 10- lo" lo' 

t [MCS] t [MCS] 



FIG. 7: (a) Number of adsorbed segments, Nadsit)), versus time t for regular A_B-copolymers with block size M = 1 64 and 
length A*' — 256. For comparison, the transient of a homopolymer is shown by a solid line too. The time interval, taken by 
the initial "shoulder", is shown in the upper left inset. The lower inset displays the variation of the scaling exponent, a, for 
the time of adsorption r oc A*'" versus block length relationship, (b) The same as in (a) but for random copolymers of length 
N — 256 and different composition p = 0.25, 0.5, 0.75. For p = 1 one has the case of a homopolymer. The inset shows the 
variation of a with p. 

The more complex adsorption kinetics, shown in Fig. [7^ for regular multiblock copolymers of block size M, and in 
Fig. [7|d for random copolymers, suggests however that the power-law character of the order parameter variation with 
time is retained except for a characteristic 'shoulder' in the adsorption transients. Indeed, one should bear in mind 
that the zipping mechanism, assumed in our theoretical treatment, is by no means self-evident when the file of sticking 
A-monomers is interrupted by neutral B segments. The characteristic shoulder in the adsorption transients of regular 
multiblock copolymers manifests itself in the early stage of adsorption and lasts progressively longer when M grows. 
The temporal length of this shoulder reflects the time it takes for a segment from the second adsorptive A-block in 
the polymer chain to be captured by the attractive surface, once the first A-block has been entirely adsorbed. For 
sufficiently large blocks one would therefore expect that this time interval, ti, associated with the capture event, 
will scale as the Rouse time, M^+^'', of a non-adsorbing tethered chain of length M. The observed ri versus M 
relationship has been shown in the upper left inset in Fig. [7^. The slope of « 1.49 is less that the Rouse time scaling 
exponent, 2.18, which one may attribute to the rather small values of the block length M that were accessible in our 
simulation. One should also allow for scatter in the end time of the shoulder due to the mismatch in the capture 
times of all the successive A-blocks in the course of our statistical everaging over many chains during the computer 
experiment. 

Somewhat surprisingly, a which describes the scaling of the total adsorption time with polymer size, r cx iV", is 
observed to decline as the block size M is increased - in contrast to the general trend of regular multiblock copolymers 
which resemble more and more homopolymers (where a ^ 1 + ly), as the block size M — > oo. Evidently, the frequent 
disruption of the zipping process for smaller blocks M slows down the overall adsorption. 

In the case of random copolymers. Fig. [TJs, the transients resemble largely those of a homopolymer chain with the 
same number of beads again, apart from the expected difference in the plateau height which is determined by the 
equilibrium number of adsorbed monomers. A rescaling of the vertical axis with the fraction of sticking monomers, 
p, however, does not lead to coinciding plateau heights - evidently the loops, whose size also depends on p, affect the 
equilibrium number of adsorbed monomers. The variation of the observed scaling exponent a with composition p is 
shown in the inset to Fig. [7|d wherefrom one gets a « 1.6 with a being largely independent of p. Note that this value 
is considerably lower than the power of 2.24 which has been observed earlier [2^, however, for very short chains with 
only 10 sticking beads. One may conclude that even for random copolymer adsorption the typical time of the process 
scales as r oc TV", as observed for homo - and regular block copolymers. It is conceivable, therefore, that an effective 
zipping mechanism in terms of renormalized segments, that is, segments consisting of an A and B diblock unit of 
length 2M for regular multiblock copolymers provides an adequate notion of the way the adsorption kinetics may be 
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treated even in such more complicated cases. For random copolymers the role of the block length M would then be 
played by the typical correlation length. 



D. Time evolution of the distribution functions - MC data 



One gains most comprehensive information regarding the adsorption process from the time evolution of the different 
building blocks (trains, loops, and tails) Probability Distribution Functions (PDF) 39]. From the MC simulation data, 
displayed in Fig. [5^, for example, one may verify that the resulting distribution D{h,t) of different train lengths is 
found to be exponential, in close agreement with the theoretically expected shape fs^, predicted under the assumption 
that local equilibrium of loops of unit length is established much faster than the characteristic time of adsorption itself. 
When scaled with the mean train length hav{t) = {h{t)), at time t, in both cases for e/ksT — 3.0 and 5.0 one finds 
an almost perfect straight line in semi-log coordinates. One may thus conclude that D{h, t) preserves its exponential 
form during the course of the adsorption process, validating thus the conjecture of rapid local equilibrium. The latter 
however is somewhat violated for the case of very strong adsorption, e/kgT = 5.0, where the rather scattered data 
suggests that the process of loop equilibration is slowed down and the aforementioned time separation is deteriorated. 

The PDF of loops W{k,t) at different times after the onset of adsorption is shown in Fig. [8}3. Evidently, the 
distribution is sharply peaked at size one whereas less than the remaining 20% of the loops are of size two. Thus the 
loops can be viewed as single thermally activated defects (vacancies) comprising a desorbed single bead with both of 
its nearest neighbors still attached to the adsorption plane. As the inset in Fig. [5Jd indicates, the PDF of loops is 
also described by an exponential function. The PDFs for loops at different time collapse on a master curve, if scaled 
appropriately with the instantaneous order parameter n{t)/N. Eventually, in Fig. [9^ we present the observed PDF 
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FIG. 8: (a) Distribution of train lengths during the adsorption process of a homopolymer chain with A'^ — 256 at two strengths 
of the adsorption potential e, shown in semi-log coordinates. PDFs for different times (in units of 10^ MCS) collapse on master 
curves when rescaled by the mean train length havit). (b) Distribution of loop lengths W{k,t) for N = 256 and e/ksT = 4.0 
during ongoing polymer adsorption. In the inset the PDF is normalized by n{t) and shown to be a straight line in log-log 
coordinates. 



T{l,t) of tails for different times t after the start of adsorption, and compare the simulation results with those from 
the numeric solution of Eq. (|4.ip . taking into account that T{1, t) = P{N — l,t). One may readily verify from Fig. [9] 
that the similarity between simulational and theoretic results is really strong. In both cases one starts at t = 1 with 
a sharply peaked PDF at the full tail length l{t = 1) = N. As time proceeds, the distribution becomes broader and 
its maximum shifts to smaller values. At late times the moving peak shrinks again and the tail either vanishes, or 
reduces to a size of single segment which is expressed by the sharp peak at the origin of the abscissa. 



Summary 



The main focus of this contribution has been aimed at the adsorption transition of random and regular multi- 
block copolymers on a flat structureless substrate whereby by different means - scaling considerations and computer 
simulations - a consistent picture of the macromolecule behavior at criticality is derived. 
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FIG. 9: (a) Distribution of tail size for different times (in units of lO' MCS) during the polymer chain adsorption for a chain 
with A'^ = 256 at e/ksT — 4.0. (b) The same as in (a) as derived from the solution of the ME for chain length A'^ = 32. For 
better visibility the time slices for t — 1, 5, 30 100, 150, 200, and 300 are shifted along the time axis and arranged such that 
the initial distribution for t = 1 is represented by the most distant slice. 

As a central result one should point out the phase diagram of regular multiblock adsorption which gives the increase 
of the critical adsorption potential with decreasing length M of the adsorbing blocks. For very large block length, 
— > 0, one finds that the CAP approaches systematically that of a homogeneous polymer. 

The phase diagram for random copolymers with quenched disorder which gives the change in the critical adsorption 
potential, e^, with changing percentage of the sticking ^-monomers, p, is observed to be in perfect agreement with the 
theoretically predicted result which has been derived by treating the adsorption transition in terms of the "annealed 
disorder" approximation. 

Evidently, a consistent picture of how some basic polymer chain properties of interest such as the gyration radius 
components perpendicular and parallel to the substrate, or the fraction of adsorbed monomers at criticality, scale when 
a chain undergoes an adsorption transition emerges regardless of the particular simulation approach. An important 
conclusion thereby concerns the value of the universal crossover exponent = 0.5 which is found to remain unchanged, 
regardless of whether homo-, regular multiblock-, or random polymers are concerned. Thus the universality class of 
the adsorption transition of a heteropolymer is the same as that of a homopolymer. 

Concerning the adsorption kinetics of a single polymer chain on a flat surface, it is shown that within the "stem- 
flower" model and the assumption that the segment attachment process follows a "zipping" mechanism, one may 
adequately describe the time evolution of the adsorbed fraction of monomers and of the probability distribution 
functions of the various structural building units (trains, loops, tails) during the adsorption process. For regular 
multiblock and random copolymers it is found that the adsorption kinetics strongly resembles that of homopolymers. 
The observed deviations from the latter suggest plausible interpretations in terms of polymer dynamics, however, it 
is clear that additional investigations will be warranted before a complete picture of the adsorption kinetics in this 
case is established too. 
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